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Abstract 

1 Coastal transport in the Bay of Palma, a small region in the island of Mallorca, 

2 Spain, is characterized in terms of Lagrangian descriptors. The data sets used for 

3 this study are the output for two months of a high resolution realistic numerical 

4 model, ROMS, forced atmospherically and with a spatial resolution of 300 m. Semi- 

5 persistent Lagrangian coherent structures (LCS) are found during both months. We 

6 study the fluid interchange between the Bay and the open ocean by computing par- 

7 tide trajectories and residence times. We found clear differences between autumn 

8 and summer. Stirring is smaller in autumn (October) than in summer (July), pre- 

9 senting diurnal variability at both months, probably induced by the effect of the 

10 wind. The escape rate of particles out of the Bay is qualitatively different, with 

11 more particles escaping to the ocean in October than in July. We show that LCSs 

12 separate regions with different transport properties, and display spatial distribu- 

13 tions of residence times on synoptic Lagrangian maps. Finally, we compare the 

14 transport at the bottom layer, strongly influenced by the bottom drag, with surface 

15 transport strongly forced by the wind. 

16 1 Introduction 

17 The study of transport and mixing in coastal flows is of major interest because of their 

18 economic and ecological importance. Due to the particularities that they present, like 

19 influence of complex topography, coastline shape and the direct driving at the surface by 

20 highly variable wind forcing, coastal flow dynamics remains still poorly understood. 

21 Recently, coastal observations and modeling efforts in different regions have been ad- 

22 dressed from the Lagrangian point of view: [25j showed that Lagrangian Coherent Struc- 

23 tures (LCSs) computed from velocity fields obtained from HF Radar measurements can 

24 be used to predict pollutant dispersion in the coast of Florida; [13] detected LCSs with 

25 HF Radar data in the Gulf of Eliat, Israel, and [18] studied small-scale properties of 

26 dispersion measurements obtained from HF Radar data in the Gulf of La Spezia, Italy. 

27 Besides Radar measurements, LCSs obtained from velocity data of high resolution nu- 
2B merical models have been used to analyze the effect of the waves on LCS in the Bay of 
29 Palma, Spain [12], to study the transport in the tidal flow of Ria de Vigo, Spain [20], or 



30 to study the water quality of a very small coastal region, the Hobie Beach, USA [TT] Also, 

31 data from drifters released in the Santa Barbara Channel were used by [30] to characterize 

32 relative dispersion, very useful to improve Lagrangian stochastic models. The application 

33 of Lagrangian techniques to study the dynamics in a shallow lake (small closed basin) has 

34 been performed in j32j . 

35 In this work we study some transport properties in the Bay of Palma using Lagrangian 

36 techniques. The Bay is a semi-enclosed basin located in the southwest of the island of 

37 Mallorca (western Mediterranean sea), whose coastal flow is strongly induced by the wind, 

38 particularly characterized by the presence of persistent sea breeze during the summer 

39 months. Forcing by tides is almost negligible with a tidal amplitude of less than 0.25r7i. 

40 The Lagrangian diagnosis will be obtained from velocity data of a realistic numerical 

41 model at high resolution, which resolves spatial scales of a few hundred of meters. We 

42 investigate the horizontal transport during two months corresponding to different seasons 

43 (autumn and summer) at the sea surface and at the deepest layer, in order to highlight 

44 the effect of the wind on transport. We compute the barriers and avenues to transport 

45 (LCS) from lines of high values of Finite-Size Lyapunov Exponents (FSLE). We also 

46 present calculations of residence times and show synoptic Lagrangian maps (SLM) of these 

47 times [26] , which will allow us a detailed visualization of the interchange of fluid particles 
4B between the Bay and the open sea. As a central point in our study, the relationship 
49 between LCSs and areas of different residence times will be analyzed. 

so The organization of this paper is as follows. The data set used in the computations 

51 and the area of study is described in Section [2j Section [3] presents a brief overview of 

52 the Lagrangian tools that are used. Before presenting the Lagrangian results, we show in 

53 Section |4] a short list of Eulerian results by studying the velocties in the Bay. We show 

54 in Section [5] a coastal characterization of stirring of the Bay of Palma in terms of FSLE 



55 and residence times. Using the definition of LCS given in Section |3} Lagrangian barriers 

56 are identified in the domain of interest. We compute escape rates and residence times 

67 of fluid particles to describe the transport relation between the Bay and the open ocean. 

68 We provide possible mechanisms to explain differences in the residences times and FSLE 
59 between different seasonal months. Finally we summarize the main results in Section [6j 

60 2 Data and characteristics of the study region 

61 2.1 Area of study 

62 The island of Mallorca (Fig. fl]) is part of the Balearic Islands Archipelago and is located 

63 in the center of the western Mediterranean (between 39° and 40°N and 2.50° and 3.50° 

64 E). The Bay of Palma is a nearly semi-circular and semi-enclosed basin located in the 

65 southwest coast of Mallorca and it can reach depths of more than 60 m. The Bay of 

66 Palma is defined as the water mass inside the square in Fig|T| consisting of a northern 

67 limit at 39°34'N, a southern limit at 39°24'N, and 2°30'E and 2°45'E as the western and 

68 eastern limits, respectively. The open boundary to the sea is in the southern part and it 

69 is 20 km wide. 

70 The size of the Bay is smaller than the Rossby radius of deformation at these latitudes, 

71 and the main circulation is determined by the bathjTiietry at the bottom layer and by the 

72 wind and sea breeze at the surface layer. During summer there are persistent sea breeze 

73 conditions. In July and August, the weather is often almost identical from one day to the 

74 next. In the vicinity of the Bay and along the southern coast of Mallorca the breeze blows 

75 from the south-west. Several studies [331 El], have pointed out that the meteorological 

76 conditions of Mallorca (intense solar radiation, clear skies, soil water deficit, dryness, weak 

77 surface pressure gradients, etc.) favors the development of sea breeze, often from April to 

78 October, and almost every day during July and August. 



79 2.2 Data 

80 The velocity data sets were obtained from the numerical model ROMS (Regional Ocean 

81 Model System). ROMS is a free surface, hydrostatic, primitive equation ocean model. 

82 The model uses a stretched, generalized nonlinear coordinate system to follow bottom 

83 topography in the vertical, and orthogonal curvilinear coordinates in the horizontal [361 

84 [13] . At each grid point, horizontal resolution Aq is the same in both the longitudinal, 0, 

85 and latitudinal, 6*, directions. 

86 We run the simulation with a resolution of Aq = 0.0027° (~300r77,, ROMS300), which 

87 is itself nested into a larger and coarser grid with Aq =1/74° (~1500r7i). The ROMS300 

88 domain covers 39°12'N - 39°36'N (latitude), and 2°24'E - 3°6'E (longitude). The total 

89 number of grid nodes is 260 x 148. Vertical resolution is variable with 10 layers in total. 

90 All domains were forced using realistic winds provided by the PSU/NCAR mesoscale 

91 model MM5. The initial vertical structure of temperature and salinity was obtained from 

92 the Levitus database ^7\ d] . 

93 We will manage velocity data from the surface layer and the bottom layer for the 

94 grid of Ao ~ 300m. This domain allows us to analyze the fluid interchange between the 

95 Bay and the open ocean, using a high resolution velocity field. Vertical displacements 

96 between layers are neglected, which is justified by the small integration times we will use. 

97 Nevertheless, close to the coast they can have an impact that will be the subject of future 

98 work. We will study two different intervals of time: one starting on October 5th, 2008 

99 and finishing on October 29th, 2008; and the other extends from July 1st, 2009 until 

100 July 26th, 2009. Temporal resolution is 15 minutes and 10 minutes for October and July, 

101 respectively, resulting in a total of 2375 snapshots of the velocity field for October, and 

102 3744 for July. 



103 3 Methodology 

104 3.1 LCSs and particle dispersion from FSLE 

105 Our methodology is based on the Lagrangian analysis of coastal flows. In the Lagrangian 

106 view, particles are advected by the flow and their horizontal motion (neglecting motions 

107 between model layers) is governed by the differential equations 

— = v,^{x,y,t), (1) 

— = Vy{x,y,t), (2) 

108 where {x{t),y{t)) are the west-east and the south-north coordinates of the trajectories 

109 and {vx, Vy) are the eastwards and northwards components of the velocity. Because of the 

110 small sizes involved, we will use a Cartesian coordinate system. 

111 LCSs [ISl El ES], are roughly defined as the material lines organizing the transport 

112 in the flow. They are the analogs, for time-dependent flows, of the unstable and stable 

113 manifolds of hyperbolic fixed points. Among other approaches [2S1 123 IIZ|, ridges of the 

114 local Lyapunov Exponents provide a convenient tool to locate them. In our case, we use 

115 the so-called Finite-Size Lyapunov Exponents (FSLEs) which are the adaptation of the 

116 asymptotic classical Lyapunov Exponent to finite spatial scales 0, E]- Thus, they are 

117 specially suited for closed or semi-closed basins. FSLEs are a local measure of particle 

118 dispersion and thus of stirring and mixing, as a function of the spatial resolution, serving 

119 to isolate the different regimes corresponding to different length scales of the oceanic flows, 

120 very useful in coastal systems [7j. 

For any pair of fluid particles, one of them located at x, the FSLE at time to and at 
the spatial point x is given by the formula: 

A(x,to,5o,<5/) = ^ln/, (3) 

\t\ Oq 
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121 where So is the initial distance of the two given particles, and Sf is their final distance. 

122 Thus, to compute the FSLEs we need to calculate the time, r, needed for the two particles 

123 initially separated 6q, to get a final distance 6f (in this way the FSLE represents the inverse 

124 time scale for mixing up fiuid parcels between length scales 6o and 6f). To obtain this 

125 time we need to know the trajectories of the particles (from Eqs. ([I]) and (|2|) which gives 

126 the Lagrangian character to this quantity. The FSLEs are computed for the points x of 

127 a lattice with lattice spacing coincident with the initial separation of fiuid particles ^o- 

128 Numerically we integrate the equations of motion using a standard, fourth-order Runge- 

129 Kutta scheme, with an integration time step of dt = 15 minutes in October and dt = 10 

130 minutes in July. Since velocity information is provided just in a discrete space-time grid, 

131 spatiotemporal interpolation of the velocity data is achieved by bilinear interpolation. For 

132 the spatial scales that define FSLEs, we take 6f = 0.1°, i.e., final separations of about 

133 10 km, because of the size of the Bay. On the other side, we take 6^ equals to 75 m, 

134 four times smaller than the resolution of the velocity field, Aq = 300 m. Since we are 

135 interested only in fast time scales, our integrations are restricted to 5 days. Locations 

136 for which the final separation at the end of this period has not reached the prescribed 

137 6f = 10 km (or for which particles have been trapped by land) are assigned a value A = 0. 

138 FSLEs can be computed from trajectory integration backwards and forwards in time. 

139 Their highest values as a function of the initial location, x, organize in filamental struc- 

140 tures approximating relevant manifolds: ridges in the spatial distribution of backward 

141 (forward) FSLEs identify regions of locally maximum compression (separation), approx- 

142 imating attracting (repelling) material lines or unstable (stable) manifolds of hyperbolic 

143 trajectories, which can be identified with the LCSs [IHl El EHl [3HI [19], and characterize 

144 the fiow from the Lagrangian point of view [211 [22] • Attracting LCSs associated to back- 
us ward integration (the unstable manifolds) have a direct physical interpretation [2T1 [8] |9]. 



146 Tracers (chlorophyll, temperature, ...) spread along these attracting LCSs, thus creating 

147 their typical filamental structure [STJ EH EU E] . 

148 We close this section by noting that the relationship between LCSs and Lyapunov 

149 exponents is based on heuristic arguments which may not be correct in some cases (see 

150 for example [IS])- The comparisons between FSLE and residence time maps presented 

151 later support the validity of the FSLE approach in the present situation. 

152 3.2 Escape and residence times 

153 Another characteristic time-scale for transport processes in open flows is the so-called 

154 escape rate [23]. This quantity measures how quickly particles trajectories escape from a 

155 domain. If we initiate A^(0) particles in a flow, we can measure how the trajectories escape 

156 the pre-selected region. In the case in which the decay in the number of particles remaining 

157 in the region up to time t, N{t), decays exponentially with time, N{t)/N{0) ~ e"''*, there 
15B is a well defined escape time defined as the inverse of the escape rate k: t^ = I/k. For 

159 the range of times explored in our work, we will see that the particle escape is close to 

160 exponential and then we can estimate the value of Xg. 

161 Te is a global quantity associated to the whole basin. A more detailed description of 

162 the transport processes can be obtained by other suitable Lagrangian quantities such as 

163 residence times [H El [ini El]. The particle residence time (RT) is defined as the interval 

164 of time that a fluid particle remains in a region before crossing a particular boundary. 

165 For each fluid particle inside the Bay at an initial time, we need to compute two times: 

166 the forward exit time, tf, computed as the time needed for a particle to cross the line 

167 delimiting the Bay, taking the forward-in-time dynamics; and the backward exit time, 

168 tb, the same but in the backward-in-time dynamics. The residence time is defined as 

169 RT = tf + th- RTs can be displayed in plots named Lagrangian Synoptic Maps [2S], in 
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170 which the residence time of each fluid particle is referenced to its initial position on the 

171 grid. 

172 4 Preliminary Eulerian description 

173 A first approach to the transport process in the Bay can be a description from the Eulerian 

174 point of view, by studying averages of the velocity field. To do this we consider separately 

175 the meridional Vy and zonal v^ components of the surface flow, and we analyze the time 

176 evolution of their spatial averages. 

177 Figures |2] a) and b) show the time series of data taken every 6 hours (solid line), and 

178 daily average time series (dashed-line) of Vx and Vy for October and July, respectively. 

179 During the two months both components of the flow present daily variability related to the 

180 presence of land and sea breezes. In July the meridional fluctuations are more noticeable, 

181 being (vy) very small. In October there are some large fluctuations of low frequency, 

182 probably induced when the large-scale winds are variable. 

183 Comparing the velocity components of both months we observe quantitative differ- 

184 ences. The values of Vy in the case of October range from -0.6 to 0.6 m/s (bottom panel 

185 of Fig. |2]a), while in the case of July, Vy is two orders of magnitude smaller, ranging from 

186 -0.004 to 0.008 m/s (bottom panel of Fig. |2]b). On other hand, Vx are similar during 

187 October and July. In October, Vx ranges from -0.6 to 0.3 m/s (top panel of Fig. [2] a), 

188 the same order of magnitude than the meridional velocity, resulting in circular motions 

189 (clockwise along the Bay). In July the situation is signiflcantly different. The zonal veloc- 

190 ity ranges from -0.8 to 0.2 m/s (top panel of Fig. b), much larger than the meridional 

191 velocity, resulting in a flow dominated by the zonal component. 



192 5 Lagrangian Results 

193 5.1 Average characterization of stirring 

194 We now describe the Lagrangian results. First we compute the temporal average (over 

195 the months of October and July) of the FSLEs for the surface layer, and for July in 

196 the bottom layer. This calculation helps us to unveil areas of different stirring and the 

197 differences between layers and months. 

19B The surface computations for the different seasonal months, October and July (Fig. Is] 

199 a, b) show different values and spatial distributions of stirring. The Bay of Palma appears 

200 to be an area with important activity. Average FSLE field looks more homogeneous in 

201 July than in October. During October filamental structures of high values of FSLE are 

202 accumulated over the northeast side of the Bay, forming a linear structure running from 

203 north to south-east which can act as a barrier, and therefore dividing the Bay in two 

204 flow regions of qualitatively different dynamics. The response of the surface currents to 

205 forcing by the wind could be an explanation of the difference between seasons. The effect 

206 of the terrain topography on stirring is clear in Fig. Islc), where FSLEs are computed at 

207 the deepest layer for July. The high values of time-averaged FSLEs are located close to 

208 a region of high bathymetry gradient. This feature is likely owing to the presence of a 

209 region of compression of the flow induced by the bathymetry. 

210 5.2 Coastal LCSs 



211 The temporal averages computed in Sec. 5.1 give us a rough idea of stirring in the Bay. 



212 More detailed information is obtained by looking at non-averaged quantities, that may 

213 reveal the existence of barriers to transport. Figure |4] shows the location of the high 

214 backward FSLE values (LCSs), appearing as a network of lines, computed at successive 

215 instants of time in October. These temporary structures can remain for one or more days, 

10 



216 as happens in October, or they can appear in the same location periodically (not shown). 

217 We stress here the appearance of a clear barrier, from north to south-east, that divides 

218 the Bay in two areas that correlates with the temporal average in Fig. [3] a). This barrier 

219 appears in almost the same position in different days, remaining without displacing too 

220 much. To effectively see that it acts as a barrier we have considered the evolution of 

221 virtual particles released at both sides of the barrier. Red and black particles do not mix 

222 and they tend to spread along the barrier (confirming that, as expected, it is an attracting 

223 line) . 

224 5.3 Transport between the Bay of Palma and the open sea 

225 In this section we study the surface transport of particles in and out of the Bay. To have an 

226 idea of the time scales involved in this interchange we proceed by computing the number 

227 of particles remaining in the Bay, N{t), averaged over different starting times (separated 
22B by IShours in order to collect the information of diurnal and nocturnal signal; this gives 

229 us 15 different simulations to be averaged for each month) as a function of the integration 

230 time t. A particle is considered to leave the Bay when crossing the red open-sea boundary 

231 in Fig. [T| so that particles landing on the coast are considered not escaped. Fig. |5] shows 

232 the different average decays for October and July. In both cases N{t) is reasonably fitted 

233 by an exponential in the considered time-range, thus identifying the escape rates k = 0.62 

234 and 0.47 days~^, respectively. The corresponding escape times, given by the inverse of 

235 the escape rate, are, respectively, t^. = 1.61 and 2.12 days. The exchange of fluid particles 

236 between the open ocean and the Bay is then less active in summer than in autumn. 

237 Next we compute synoptic maps of the residence times. As was indicated in Sec. 



3.2 the residence time of the particles throughout the study area is considered as the 



239 sum of the entry time (tb) and the escape time (t/). To compute tf and tb particles are 
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240 initialized every 6 hours in a regular grid of 75m spacing and they are integrated forward 

241 and backward in time during 5 days. We consider that 5 days is a proper integration time 

242 according with the time scales associated with the coastal processes of this small Bay, 

243 and also owing to the short period of the available data. In these computations we assign 

244 the maximum possible value of tf and tf, (5 days) to the fluid particles that remain in the 

245 pre-selected area after the 5 days of integration. 

246 In Fig. [6] we color the initial positions of particles in the Bay attending to the time 

247 they transit through the Bay {RT = tf + tf,) for different days. Initial positions of particles 
24B with short residence times are indicated in blue in Fig. [6] Regions from where particles 

249 have longer residence time (i.e. take more time between entry and escape) are marked in 

250 red/brown. 

251 These maps show that the spatial distribution of residence times of particles can be 

252 complex and time depending, presenting different patterns at different times. A number 

253 of small structures can be observed, including thin filaments or small lobes. Comparing 

254 both months, one can see differences in the RT distributions. The most noticeable is the 

255 approximate east-west alignment of the zones of similar RT in July, which is not seen 

256 in October. Also, in October the values of residence times of the particles are smaller, 

257 in agreement with the global rate estimations showed before. A common feature is the 
25B southwest region with low values of residence times, because in this region there are not 

259 coastal boundaries and it is totally open to the ocean. 

260 In order to reveal regions with different persistent transport properties we compute 

261 time averages of the spatial distributions of residence times. We average 6-hour ly snap- 

262 shots of RT during 15 days (i.e. 60 snapshots) for each month. The results, plotted in 

263 Fig. [TJa) and b), show the common features that, in general, the low values of RT are for 

264 particles initiated close to the open ocean, specially in the southwest part, and high values 

12 



265 are for particles started near the coast, as expected. However, on average, the residence 

266 time is larger in July than in October (3.25 days in July and 1.51 days in October) con- 

267 sistent with the behavior of the corresponding values of Xg. Also, in July there is a clear 
26B boundary between the interior of the Bay to the north, with large average residence times 

269 and the open sea to the south, whereas the boundary between high and low residence 

270 times in October is well inside the Bay, aligned with the Lagrangian structure identified 

271 from the FSLE analysis, as will be discussed in the next section. 

272 Another feature observed in movies of particle trajectories (not shown) is that in 

273 October fluid particles tend to circulate mostly clockwise, while in July they are oscillating 

274 along the zonal direction (see Section ffl). This difference is likely to be responsible for the 

275 different behavior between both months. 

276 5.4 Relation between LCSs and RTs 

277 We now examine the connection between regions of different residence times with LCSs. 

278 To compare RT and FSLE we have superimposed in Fig. [6] the fllaments of high values 

279 of forward (white) and backward FSLE (black) values on the spatial distribution of res- 

280 idence times. This Figure shows a clear correspondence between structures of RT and 

281 FSLE. The LCSs given by high values of FSLE clearly separate regions with different 

282 values of residence times, conflrming the value of the FSLE technique to identify bound- 

283 aries between different flow regions and barriers to transport. [32j studied the relation 

284 between residence time and FSLE for a wind-forced hydrodynamical model of a shallow 

285 lake. They found that areas with long residence time visualize the stable manifolds of 

286 the so-called chaotic saddle, a structure controlling the escape properties at long times. 

287 In our case, our integrations are restricted to times too short to characterize long-time 

288 chaotic behavior, but still there is a good correspondence between the FSLE structures 
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289 characterizing divergent or convergent trajectories, and escape or residence times. 

290 Fig. [7] a) shows a time average over October of the spatial distribution of RT, to be 

291 compared with the corresponding average figure (Fig. [3] a) for FSLE. It is evident that the 

292 region in the north and east side of the Bay with high values of RT is separated from the 

293 rest by a region of high values of FSLE. This can be explained by the presence of persistent 

294 barriers which do not allow particles to escape from the northeast side of the Bay, and 

295 thus separating the Bay in regions with different residence times. In July the situation is 

296 different, because the spatial distribution of FSLE (Figjsjb) and RT (Fig. [7|b) is almost 

297 homogeneous, with higher values over the whole area of the Bay, and lower values in the 
29B small region bordering the open ocean. Nevertheless, the instantaneous configurations 

299 of high FSLE lines (Fig. IgI) reveal their general east-west alignment, consistent with the 

300 orientation of the boundaries between areas of different RTs. 

301 5.5 Variability of RT and FSLE 

302 The differences of residence times and FSLEs in the two considered months indicate that 

303 the dynamics of the flow is qualitatively different. Now we analyze the time evolution of 

304 their spatial averages. 

305 Figures [s] a) and b) show the time series of the spatial mean of residence times (top 

306 panel) and FSLE (bottom panel) for October and July, respectively. 

307 Comparison between time evolution of spatial average of the RT for the different 

308 months confirms, again, that particles tend to stay longer times in the Bay during July 

309 than in October. The values of RT vary approximately from 0.25 days to 3 days in October 

310 (Fig. |8]a, top), and from 2 days to 6 days in July (Fig. |8]b, top). The same happens with 

311 FSLE, higher values correspond to July and lower ones to October. Diurnal fluctuations, 

312 likely related to the effect of the sea breeze, are evident in RT and FSLE for both months. 
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313 In October there are some large fluctuations of low frequency in RT, probably induced by 

314 the variability of remote forcing winds. On the other hand, during October, minima of 

315 RT correspond to maxima of FSLE, and maxima of RT correspond to minima of FSLE. 

316 In July, the relationship between FSLE and RT is looser and only observed in the high- 

317 frequency fluctuations. This difference in behavior probably arises from the fact that high 
31B values of FSLE determine clear and well-defined barriers to transport only in the case of 



319 October (see Sect. 5.3 and 5.4). In July lines of high FSLE remain nearly horizontal and 

320 probably represent regions of high zonal shear everywhere in the Bay, that by themselves 

321 do not imply stronger or weaker escape towards the South. 

322 5.6 Transport at the bottom layer 

323 In this subsection we compare the main Lagrangian characteristics at the bottom layer, not 

324 driven directly by wind, with those at the surface. In Fig. [91a) we show an instantaneous 

325 map of the residence times in the bottom layer one day of July, overlayed with lines of 

326 high FSLE values. Again, the spatial distribution is inhomogeneous, and we find high 

327 values of RT over all the Bay except very close to the ocean. The correlation of RT values 

328 with FSLE lines is weaker than in the upper layer, but still we see that the relatively 

329 lower values of RT in the western part of the Bay at that particular day appear bounded 

330 by backwards FSLE lines, indicating a temporal escape route of particles in that region 

331 towards the southwest. The spatial distribution of the time average of RT plotted in 

332 Fig. M b) shows that highest average values of RT are concentrated in the northwestern 

333 region of the Bay. Fig. [s] c) displays high values of FSLE located precisely in the same 

334 region where the RT qualitatively change to high values. This suggests the presence of 

335 persistent barriers that separate this southeastern region from the rest in this bottom 

336 layer. The formation of these persistent LCSs is, probably, produced by the gradient of 
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337 the bathymetry (see Fig. [I]). 

33B The time evolution of the spatial average of RT and FSLE are plotted in top and 



339 bottom panel of Fig. 10 respectively. Contrarily to the surface, in the bottom layer the 

340 diurnal fluctuations in the time series of RT disappear, showing that the flow at this depth 

341 is not directly influenced by breeze. The RT values are larger than in the surface, and 

342 therefore the interchange between the ocean and the Bay is less intense at bottom layers. 

343 This is a consequence of the smoothness of the flow produced by the absence of the wind 

344 forcing at deepest levels. In average, there is a negative relation between stirring and 

345 residence times: when the flow is more dispersive the particles transit during less time 

346 over the Bay, so that maxima of RT correspond to minima of FSLE and vice- versa. 

347 6 Conclusions. 

34B Properties of coastal transport in the Bay of Palma, which is a small semi-enclosed region 

349 of the Island of Mallorca, were studied in a Lagrangian framework, by using a realistic 

350 model velocity data at high resolution. We have applied two Lagrangian methods (FSLEs 

351 and RT) to analyze the small scales of these coastal currents. LCSs have been detected as 
362 ridges of FSLE, and virtual experiments with particle trajectories have shown that these 

353 structures really act as barriers, organizing the coastal flow. Global and average aspects 

354 of the transport in different seasonal months show that, in the period studied, in autumn 

355 there is more interaction between the Bay and the open ocean than in summer. One 

356 explanation is that during July the flow conditions restrict motion of the coastal marine 

357 surface to the zonal direction, preventing the flow to enter or scape toward the open 
35B ocean. The transport of particles at the deepest layer, strongly influenced by the bottom 

359 topography, is less active than at the surface. Regions with different values of RT are 

360 separated by ridges of FSLE, proving the fact that FSLE separate regions of qualitatively 
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361 diflFerent dynamics also in small coastal regions. Thus, we think that these Lagrangian 

362 quantities can be used as key variables determining the dynamics and health of other 

363 bays or estuaries, particularly in relation with human activities. Future improvements 

364 include the adaptation of these methods to three-dimensional spaces and capture three- 

365 dimensional effects, such as upwelling and downwelling in coastal areas, and analyzing 

366 longer periods of time. 
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Figure 1: Bathymetry contours (in meters) of the model domain. The black box indicates 
the Palma Bay and the inset graphics give the geographical location of Mallorca Island 
in the western Mediterranean Sea. 
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Figure 2: a) Time series throughout October of the zonal (top) and meridional (bottom) 
components of the velocity field every 6 hours (solid line) and daily averaged (dashed-line); 
b) the same that a) but for July. 
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Figure 3: Spatial distribution of the time average of 6-hourly FSLEs maps over different 
months and at different layers: a) October at surface layer, b) July at surface layer, c) 
July at bottom layer. 
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Figure 4: Evolution of the locations of two sets of particles in the Palma Bay during a night 
of October, superimposed on the spatial distributions of high values of backward FSLEs. 
The colorbars specify FSLE values in days~^. Zero FSLE values, displayed as white, are 
assigned to locations for which the particles do not attain the prescribed Sf = 10km 
separation after 5-days integration. Note the highest values of FSLEs (green lines) act as 
a barrier practically dividing the Bay in two parts. The two sets of particles are deployed 
from both sides of the barrier, (a) Initial conditions of the particles on October 8 at 20:00 
GMT, 2008.; (b) October 9 at 00:00 GMT, 2008 ; (c) October 9 at 04:00 GMT, 2008 ; 
(d) October 9 at 08:00 GMT, 2008. Particles marked by black dots were released in the 
right side (northeast) of the barrier while the particles marked with red were released on 
the left side of the barrier. 
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Figure 5: Average of 15 subsequent (started at to values separated 18 hours) estimations 
of N[t), the number of particles remaining in Palma Bay at least for a lapse of time t after 
release at to- Black and red lines are for surface layer in October and July respectively. 
Dashed lines are the measured averages, and the solid lines are the indicated exponential 
fits. 
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Figure 6: Lines are the locations of top values of FSLE (greater than 0.5days~^ in October 
and greater than 1.5days~^ in July). Backward FSLE lines are colored in black and 
forward FSLEs in white. They are superimposed on spatial distributions of residence 
times in Palma Bay for different dates. The colorbars give the residence times in days, 
a) and b) correspond to two different days in October, and c) and d) in July. 
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Figure 7: Spatial distributions of time averages of 60 snapshots of 6-hourly RT values 
collected over 15 days in Palnia Bay for a) October and b) July. The colorbar units are 
days. 
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Figure 8: a) Time series throughout October of the spatial average of residence times 
(top) and spatial average of FSLE (bottom) for the surface layer of Palma Bay. b) the 
same that a) but for July. 
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Figure 9: a) Snapshot of spatial distributions of residence times at the bottom layer in 
Palma Bay corresponding to July 17, 2009 at 18:00h (GMT). Lines are the locations of 
top values of FSLE (greater than O.Sdays^^). Backward FSLE lines are colored in black 
and forward FSLEs in white, b) spatial distribution of time average of 60 snapshots of 
6-hourly maps of RT collected over 15 days in July at the bottom layer. 
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Figure 10: Time series throughout July of (top) spatial average of residence times and 
(bottom) spatial average of FSLE computed for the bottom layer. 
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